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In a recent publication [Ferreira et al, Pliys. Rev. D89 (2014) 083011] we tested the consistency of 
current astrophysical tests of the stability of the fine-structure constant a and the proton-to-electron 
mass ratio pi = m p /m e (mostly obtained in the optical/ultraviolet) with combined measurements 
of a, p and the proton gyromagnetic ratio g p (mostly in the radio band). Given the significant 
observational progress made in the past year, we now revisit and update this analysis. We find 
that apparent inconsistencies, at about the two-sigma level, persist and are in some cases enhanced, 
especially for matter era measurements (corresponding to redshifts z > 1). Although hidden sys- 
tematics may be the more plausible explanation, we briefly highlight the importance of clarifying 
this issue, which is within the reach of state-of-the art observational facilities such as ALMA and 
ESPRESSO. 

PACS numbers: 04.50.Kd, 98.80.-k 


I. INTRODUCTION 

Tests of the stability of nature’s fundamental couplings 
are becoming an increasingly important probe of physics 
beyond the standard model [1:, i2| ■ Very tight constraints 
stem from local laboratory tests using atomic clocks @, 
while astrophysical measurements allow a large lever arm 
which can probe the dynamics of the new degrees of free¬ 
dom responsible for such putative variations, and poten¬ 
tially shed light on the dark energy enigma [I]. There 
have been recent indications of possible variations Q, 
which a dedicated Large Program at ESO’s Very Large 
Telescope is aiming to test [6H8|. 

Direct astrophysical measurements of the fine- 
structure constant a and the proton-to-electron mass ra¬ 
tio pi = m p /m e are, in most cases, carried out the opti¬ 
cal/ultraviolet (there are a few exceptions to this for the 
pi case), and up to redshifts now exceeding z = 4. On 
the other hand, in the radio band, and typically at lower 
redshifts, one can measure various combinations of a, pi 
and the proton gyromagnetic ratio g p . In a recent work 
[Oj we carried out a joint statistical analysis of all existing 
data, in the context of a broad class of unification sce¬ 
narios 53. and highlighted some apparent inconsisten¬ 
cies which could be an indication that systematics may 
be affecting some of the data. 

Given the significant observational progress made in 
the past fifteen months, both in measurements of a BEl 
and in those of pi fl2Hl4| |. in this work we will update 
our previous analysis. We start in Sect. II by briefly 
summarizing the theoretical assumptions underlying the 
class of unification scenarios we will assume, and then 
provide an up-to-date list of available measurements in 
Sect. III. Our consistency analysis is then presented in 
Sect. IV, and a brief outlook follows in Sect. V. 


II. VARYING COUPLINGS AND UNIFICATION 

The Standard Model (SM) of particle physics is one of 
the most respected buildings of modern physics. The de¬ 
sire to push it forward and understand its correctness at 
sufficiently high energy scales (the Planck mass, for in¬ 
stance) raises some open issues. One hint to the fact that 
the SM might be incomplete is the running of coupling 
constants with energy. These couplings are expected to 
meet at an energy scale of Mqut ~ 10 16 GeV. Consid¬ 
ering this fact, it seems reasonable to think of a scenario 
in which all interactions stem from only one theory: a 
Grand Unified Theory (GUT). 

In a general GUT, the symmetries of the SM—the 
gauge group 51/(3) x SU( 2) x U( 1)—are unified into a 
larger symmetry group with only one coupling param¬ 
eter at the GUT energy scale Mqut ~ 10 16 GeV. This 
coupling is assumed to be an independent parameter and 
it is allowed to vary. Assuming one such GUT, one can 
relate all couplings of SM interactions to only one, fun¬ 
damental coupling. Since the latter can vary, all other 
related, lower energy, couplings must also vary and those 
variations should be related. 

There are many models that attempt to unify the SM 
interactions. We will be working with a specific class 
of models that makes the following assumptions: the 
electroweak scale is derived by dimensional transmuta¬ 
tion, all Yukawa couplings vary in the same way and the 
variation of the couplings is assumed to be a result of a 
dilaton-type scalar field. Under these assumptions it is 
possible to encapsulate the details of these class of the¬ 
ories in two parameters: R and S pim. They are 
defined as proportionality factors between variations of 
fundamental couplings and they are highly model depen¬ 
dent: 
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where R ~ 95.7 5/3A ffi. At3 and 5 = in which b t 

is the beta function coefficient for the i th SM gauge cou¬ 
pling (i = 1, 2,3) at an energy scale smaller than the uni¬ 
fication scale, A qcd is the QCD mass scale, a is the fine- 
structure constant, v is the vacuum expectation value of 
the Higgs and h stands for any of the Yukawa couplings 
(we’re assuming that all of them vary in the same way). 

Using these two parameters, one can express, in a com¬ 
pact way, the varying-coupling character of a wide class 
of theories and, by constraining them, it is possible to 
shed some light on the specificities of a putative GUT. 
Hence, in what follows, we’ll use them to relate the vari¬ 
ations of interest as [10, Ha, [13 
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^ = [0.127?- 0.05(1 + S)]—. (5) 
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Note that in this class of models the proton gyromag- 
netic ratio is less sensitive to the parameters R and S 
than the proton-to-electron mass ratio. We will use this 
observation in some of our subsequent analysis. We also 
point out that in this class of models the behavior of the 
temperature-redshift relation provides a further consis¬ 
tency test, as first discussed in 0111. 


III. CURRENT SPECTROSCOPIC 
MEASUREMENTS 

We now list the astrophysical measurements that will 
be used in our analysis. Unless otherwise stated, we will 
list them in units of parts per million (ppm). This is 
not meant to be an exhaustive list of all measurements. 
We typically use only the tightest available measurement 
for each astrophysical source. A few older measurements 
along other lines of sight have not been used, on the 
grounds that they would have no statistical weight in 
the analysis. Nevertheless, we will include some low- 
sensitivity but high-redshift measurements, as these are 
illustrative of the redshift range that may be probed by 
future facilities. As in [g], whose list we update here, 
our two exceptions regarding measurements of the same 
source are 

• Measurements using different, independent 
techniques—typically this occurs with measure¬ 
ments of /i or combined measurements using 
different molecules, and 


Object 

z 

Qab 

AQab/Qab 

Ref. 

PKS1413+135 

0.247 

a 2xl - 85 g pt i 185 

-11.8 ±4.6 

[20] 

PKS1413+135 

0.247 

a 2x1 - 57 g p g 157 

5.1 ± 12.6 

[21] 

PKS1413+135 

0.247 

a 2 g v 

-2.0 ±4.4 

[22] 

B0218+357 

0.685 

a 2 g p 

-1.6 ±5.4 

[22] 

J0134—0931 

0.765 

Q 2x 1.57^^1.57 

-5.2 ±4.3 

[23] 

J2358—1020 

1.173 

a 2 g v t\x 

1.8 ±2.7 

[24] 

J1623+0718 

1.336 

^ 9p/ 9 

-3.7 ±3.4 

M 

J2340—0053 

1.361 

a 2 g p /fi 

-1.3 ±2.0 

m 

J0501—0159 

1.561 

a 2 g p /g 

3.0 ±3.1 

m 

J1024+4709 

2.285 

a 2 fi 

100 ± 40 

[25] 

J2135—0102 

2.326 

a 2 fi 

-100 ± 100 

[25] 

J1636+6612 

2.517 

a 2 g 

-100 ± 120 

[25] 

H1413+117 

2.558 

a 2 fi 

-40 ± 80 

[25] 

J1401+0252 

2.565 

a 2 fi 

-140 ±80 

[25] 

J0911+0551 

2.796 

a 2 fi 

-6.9 ±3.7 

[26] 

J1337+3152 

3.174 

a 2 g p /fi 

—1.7 ± 1.7 

m 

APM0828+5255 

3.913 

a 2 fi 

-360 ± 90 

[25] 

MM1842+5938 

3.930 

a 2 fi 

-180 ±40 

[25] 

PSS2322+1944 

4.112 

a 2 fi 

170 ± 130 

[25] 

BR1202—0725 

4.695 

a 2 fi 

50 ± 150 

[28] 

J0918+5142 

5.245 

a 2 fi 

— 1.7 ± 8.5 

[29] 

J1148+5251 

6.420 

a 2 fi 

330 ± 250 

[28] 


TABLE I. Available measurements of several combinations of 
the dimensionless couplings a, fi and g p . Listed are, respec¬ 
tively, the object along each line of sight, the redshift of the 
measurement, the dimensionless parameter being constrained, 
the measurement itself (in parts per million), and its original 
reference. 


• Measurements obtained with different spectro¬ 
graphs . 

In these cases we do list the various available measure¬ 
ments. 

Table E contains current joint measurements of sev¬ 
eral couplings. Compared to our earlier work we have 
added the measurements of [25| , which although not par¬ 
ticularly sensitive complement the other measurements 
in terms of redshift coverage. (We thank Hugo Mes- 
sias for bringing this reference to our attention.) Note 
that for the radio source PKS1413+135 the three avail¬ 
able measurements are sufficient to yield individual con¬ 
straints on the variations of the three quantities at red¬ 
shift 2 = 0.247. This analysis was done in 0], yielding 
a null result at the two sigma confidence level. 

Table [H| contains individual a measurements. Con¬ 
servatively we only list measurements where data was 
acquired specifically for this purpose—but these are by 
now, in most cases, the ones with the smallest uncertain- 
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Object 

z 

A a/a (ppm) 

Spectrograph 

Ref. 

3 sources 

1.08 

4.3 ±3.4 

HIRES 

[UJ 

HS1549±1919 

1.14 

-7.5 ±5.5 

UVES/HIRES/HDS 

m 

HE0515—4414 

1.15 

-0.1 ±1.8 

UVES 

Sii 

HE0515—4414 

1.15 

0.5 ±2.4 

HARPS/UVES 

m 

HS1549±1919 

1.34 

-0.7 ±6.6 

UVES/HIRES/HDS 

m 

HE0001—2340 

1.58 

— 1.5 ± 2.6 

UVES 

[33] 

HE1104—1805A 

1.66 

-4.7 ±5.3 

HIRES 

UJJ 

HE2217—2818 

1.69 

1.3 ±2.6 

UVES 

[61 

HS1946±7658 

1.74 

-7.9 ±6.2 

HIRES 

ill 

HS1549±1919 

1.80 

-6.4 ±7.2 

UVES/HIRES/HDS 

[81 

Q1101—264 

1.84 

5.7 ±2.7 

UVES 



TABLE II. Available specific measurements of a. Listed are, 
respectively, the object along each line of sight, the redshift 
of the measurement, the measurement itself (in parts per 
million), the spectrograph, and the original reference. The 
recent UVES Large Program measurements are fh, §j. The 
first measurement is the weighted average from 8 absorbers 
in the redshift range 0.73 < z < 1.53 along the lines of sight 
of HE1104-1805A, HS1700+6416 and HS1946+7658, reported 
in [ll] without the values for individual systems. The UVES, 
HARPS, HIRES and HDS spectrographs are respectively in 
the VLT, ESO 3.6m, Keck and Subaru telescopes. 


Object 

z 

A/x//x 

Method 

Ref. 

B0218±357 

0.685 

0.74 ± 0.89 

nh 3 /hco+/hcn 

[34] 

B0218±357 

0.685 

-0.35 ±0.12 

NH 3 /CS/H 2 CO 

[35] 

PKS1830—211 

0.886 

0.08 ± 0.47 

NH 3 /HC 3 N 

[36] 

PKS1830—211 

0.886 

-1.2 ±4.5 

CH 3 NH 2 

[32 

PKS1830—211 

0.886 

-2.04 ± 0.74 

NH 3 

[38] 

PKS1830—211 

0.886 

-0.10 ±0.13 

CH 3 OH 

m 

J2123—005 

2.059 

8.5 ±4.2 

H 2 /HD (VLT) 

[40] 

J2123—005 

2.059 

5.6 ±6.2 

H 2 /HD (Keck) 

m 

HE0027—1836 

2.402 

-7.6 ± 10.2 

h 2 

m 

Q2348—Oil 

2.426 

-6.8 ±27.8 

h 2 

[42] 

Q0405—443 

2.597 

10.1 ±6.2 

h 2 

m 

J0643—504 

2.659 

7.4 ±6.7 

h 2 

m 

Q0528—250 

2.811 

0.3 ±3.7 

H 2 /HD 

m 

Q0347—383 

3.025 

2.1 ±6.0 

h 2 

[45] 

J1443±2724 

4.224 

-9.5 ±7.6 

h 2 

m 


TABLE III. Available measurements of fi. Listed are, respec¬ 
tively, the object along each line of sight, the redshift of the 
measurement, the measurement itself, the molecule(s) used, 
and the original reference. The recent LP measurement is [ 3 ] • 


ties. Compared to our previous analysis we have added 
the most recent Large Program result Q as well as the 
measurements of [11] . We note that the weighted mean 
of the measurements on the table is 

(—\ =0.37 ±0.94, (6) 

V a / wm 

and thus, for the first time, the one-sigma statistical un¬ 
certainty is below ppm (although the key issue is that of 
possible systematics). 

Table Mil contains individual /x measurements. Note 
that several different molecules can be used, and in the 
case of the gravitational lens PKS1830—211 several in¬ 
dependent measurements exist. Currently ammonia is 
the most common at low redshift, though others such as 
methanol, peroxide, hydronium and methanetiol have a 
greater potential in the context of facilities like ALMA 
[461 . At higher redshifts molecular hydrogen is the most 
common. 

The tightest available constraint was obtained in 
PKS1830—211, from observations of methanol transitions 
39]. A recent analysis leads to a nominally tighter bound 
14] , but the authors suggest caution in using that re¬ 
sult due to possible systematics; conservatively, we have 
therefore kept the earlier measurement. As for measure¬ 
ments along the line of sight towards J0643—504, we have 
similarly listed the more conservative one, from lilt on 
the table; for comparison [47| finds A fx = 12.7 ± 6.2 
ppm. 

The main update to these measurements is therefore 
the arrival of the first direct constraint on /x beyond z = 4 


0. Calculating the weighted mean of the low and high- 
redshift samples (z < 1 and z > 2, respectively, we find 

= —0.24 ± 0.09 (7) 

Low ,wm 



= 3.4 ±2.0, (8) 

High,wm 

in both cases this is a very mild evidence for a variation, 
though with different signs at high and low redshifts, 



IV. CONSISTENCY ANALYSIS 

We start by determining the individual values of the 
three couplings that provide the best fit to the data in Ta¬ 
ble |U The results are summarized in Fig. [TJ We obtained 
the following one-sigma constraints for the full table 

— = - 2 . 511:1 ( 9 ) 

a 

— = - 1 . 611 ? ( 10 ) 

/I 

^= 2 . 7 ™ ( 11 ) 

9p 

thus consistent with no variations. In addition to the full 
sample the figure also shows the results of the analysis 
for a few subsamples 
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FIG. 1. Best-fit parameters for the data in Q] One-dimensional confidence intervals for the relative variations of the three 
couplings (in ppm), for the full sample and the subsamples discussed in the text. 

Note that the vertical scale is different in both sets of panels. 




FIG. 2. Best-fit parameters for the data in[T] Two-dimensional 
one, two and three sigma confidence intervals in the R-S plane. 

• z < 0.25 (denoted Low in the plot, and containing 
only the measurements along the line-of-sight to 
PKS1413+135 discussed in [30]) versus z > 0.25 
(denoted High in the plot) 

• z < 1 versus z > 1 

• z < 2 versus z > 2 

We note that the likelihood profiles are usually non- 
Gaussian, both here and in the analyses in the rest of 
this section. This explains why some of the confidence 
intervals are asymmetric. 

Fig. [2] shows the two-dimensional confidence levels in 
the plane of the unification parameters R and S discussed 
in Sect. II. The one dimensional marginalized likelihoods 
for each of them are 


0 < 2 < 1 1 < 2 < 2 2 > 2 



Aa A a Aa 


FIG. 3. Two-dimensional likelihood in the a — g plane for the 
data in[U assuming no variations of g v , for the full table (top 
panel) and dividing the data into three redshift bins (bottom 
panel). In both cases the one, two and three sigma confidence 
intervals are plotted. 


R = 253 ±221, S'= 673 ±594. 


(12) 


























































































5 


Sample 

A a/a 

Ag/g 

2 < 1 

-0.9 ±4.2 

-2.7 ±9.6 

1 < 2 < 2 

Unconstrained 

Unconstrained 

2 > 2 

-2.3 ±2.3 

-2.8 ±4.7 

Full 

—1.4 ± 1.1 

-2.2 ±2.3 


TABLE IV. One-dimensional marginalized two-sigma con¬ 
straints for a and g, for the data in |H assuming no varia¬ 
tions of g p , for the full sample and various redshift bins. All 
constraints are in parts per million. 


These are fully consistent with the naive expectations 
suggesting typical values around R ~ 30 and S ~ 160 
Hot as well as with the results of our earlier work Q. 

Since in this class of models the variations of g p are, for 
typical values of the parameters R and S, smaller than 
those of /x, it’s interesting to repeat some of the above 
analysis on the assumption that g p does not vary. These 
results are shown in Fig. [3] and the ID marginalized 
constraints for the full sample and various subsamples 
are shown in Table hyi 

Note that in this case there are again deviations (at 
up to the two sigma level) from the null result, which 
wasn’t the case in the full analysis where g p was also 
allowed to vary. This change stems from the reduction 
in uncertainties afforded by the assumption of a fixed g p 
(the shift in best-fit values is comparatively small), and 
is therefore not statistically significant. In the redshift 
range 1 < z < 2 the individual values of both parameters 
are unconstrained by this dataset: only the combination 
a 2 //x is constrained, but this constraint does play a role 
in the full likelihood for this case. Finally, it’s also inter¬ 
esting to note that the degeneracy directions are different 
at high and low redshifts. 

We can also reanalyze the data of Table Q] using the 
weighed mean value of Table [TT] for a (ie, EqjGj) as a prior. 
The results of this analysis are summarized in Fig. [I] and 
Table El both for the full sample and for high/low red¬ 
shift subsamples. There is a mild (two-signra) preference 
for negative variations of both quantities, which mostly 
stems from the higher redshift data. As before, note 
that the degeneracy directions change with redshift; this 
partly explains why the preference for negative variations 
is much stronger for the full sample than the subsamples. 
On the other hand, it could also be an indication that as¬ 
suming a constant value for the parameters in this entire 
redshift range may be too simplistic. 

The complementary exercise can also be done: we can 
analyze the data of Table U using the weighed mean values 
of Table Hill for g at low and high redshifts (ie, Eqs[71[8]) as 
priors. The corresponding results are summarized in Fig. 
[5] and Table Ell We again considered both the full sample 
and high and low redshift subsamples. There is again a 
mild (two-sigma) preference for variations of both quan¬ 
tities, which again mostly steins from the higher redshift 
data. In this case the results are not as clear which is 




A p 


FIG. 4. Two-dimensional likelihood (one, two and three sigma 
confidence levels) in the g — g p plane for the data in Table 
H assuming the weighed mean value of Table HO for a (given 
by Eq[6]).In both panels the two sets of black contours show 
the likelihood for low and high redshift subsamples, and the 
shaded regions show the likelihood for the full sample. In the 
top panel the borderline between the two subsamples is at 
2 = 1, while in the bottom one it’s at z = 2. 


at least in part the result of a stronger degeneracy than 
in the previous case. One effect of the g priors (with 
their slight preference for negative variations) is to lead 
to a preference for positive values of g p (while a negative 
value is preferred for a). 

As an additional check, if we assume that g p does not 
vary and use the weighed mean of Table E] for — as a 
prior for Table U we now find, at the two-sigma confidence 
level, 


— = -2.9 ±2.2 (13) 

A 4 

compared to the more general analysis summarized in 
Table El the central value is in good agreement, while 
the uncertainty is naturally reduced. Conversely, if we 
assume that g p is fixed and use the weighted mean value 
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Sample 

A g/g 

A g P /g P 

0 < z < 1 

-0.1 ±6.8 

-2.5 ±4.5 

z > 1 

-3.1 ±5.9 

-4.2 ±6.7 

0 < z < 2 

-1.1 ±3.6 

-1.9 ±3.2 

z > 2 

-3.1 ±5.9 

-5.6 ±8.6 

Full 

-3.3 ±3.3 

-4.2 ±3.6 


TABLE V. One-dimensional marginalized two-sigma con¬ 
straints for g and g p , for the data in Tabic Q] assuming Eq. 
[6] as prior on a, for the full sample and various redshift bins. 
All constraints are in parts per million. 


Sample 

Aa/a 

A g P /g P 

0 < z < 1 

0.5 ±7.6 

-2.8 ± 16.6 

z > 1 

-2.9 ±3.0 

8.9 ±6.9 

0 < z < 2 

-0.1 ± 7.4 

-0.6 ± 15.2 

z > 2 

-2.9 ±3.0 

7.5 ±8.9 

Full 

-4.0 ±3.3 

6.9 ±7.5 


TABLE VI. One-dimensional marginalized two-sigma con¬ 
straints for a and g p , for the data in Table [T] assuming Eqs. 
m as low and high-redshift priors on g, for the full sample 
and various redshift bins. All constraints are in parts per 
million. 


of Table Hill for we find, still at the two-sigma level 

— = -1.1 ±0.8; (14) 

a 

again, this is consistent with the results of Table IVI1 
though with a significantly reduced uncertainty. Finally 
if we simultaneously use the weighted mean averages of 
Tables [TT] and [ITT] as priors for ^ and and use Table 
Uto infer g p , we obtain at the two-sigma confidence level 

^ = -2.1 ±2.8. (15) 

9p 

V. CONCLUSIONS 

In this work we have fully updated an earlier analy¬ 
sis |§| of the consistency of the various currently avail¬ 
able astrophysical tests of the stability of fundamental 
couplings. Various combinations of a , g and g p can be 
measured in the radio band, while a and g can also be 
individually measured, mostly through observations in 
the optical/ultraviolet region. Taken together, they span 
very broad redshift range—approximately from z = 0.2 
to z = 6.4—although the sensitivity of current measure¬ 
ments is also quite heterogeneous: from about 0.1 ppm 
to more than 100 ppm. 

Our results show that a joint analysis of all the com¬ 
bined measurements is consistent with no variations. 




Aa 


FIG. 5. Two-dimensional likelihood (one, two and three sigma 
confidence levels) in the a — g p plane for the data in Table [T] 
assuming the weighed mean values of Tablc lIIIl at low and high 
redshifts for g (cf. Eqs. I7l8l) .ln both panels the two sets of 
black contours show the likelihood for low and high redshift 
subsamples, and the shaded regions show the likelihood for 
the full sample. In the top panel the borderline between the 
two subsamples is at z = 1, while in the bottom one it’s at 
z = 2. 

However, joint analyses of the radio and optical/UV data 
tend to lead to mild evidence (specifically, at two-sigma 
confidence level) for variations at the few ppm level. This 
is especially the case at redshifts z > 1, corresponding to 
measurements in the matter era. 

It is noteworthy that the sensitivity with which each 
parameter can be constrained is redshift dependent. The 
best direct individual constraints on a (at the few ppm 
level) are currently in the redshift range 1 < z < 2, while 
those for g are either at z > 2 (from molecular Hydrogen 
measurements) or at z < 1 (from measurements using 
other molecules, such as ammonia or methanol). Simi¬ 
larly, different products of a, g and g p can be constrained 
at different redshift. This is relevant for optimizing ob¬ 
servational strategies and selecting targets for future fa¬ 
cilities. 

Admittedly the evidence for possible variations thus 
inferred is not strong, and the most likely explanation 
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for them is that systematic errors have not been fully 
accounted for. In the case of the optical measurements, 
possible sources for these have recently been studied in 
some detail [fH-lsl l48|. In any case, clarifying these dis¬ 
crepancies is essential. Finding lines of sight where these 
measurements can be carried out both in the optical and 
in the radio bands would provide an ideal test of possible 
systematics, but the number of such targets is likely to 
be small. 

A more immediate goal will be to extend the range of 
redshifts where sensitive measurements of a and fj, can 
be made. The imminent arrival of observational facilities 
such as ALMA and ESPRESSO will make this possible. 


This will signal the start of a new era of precision consis¬ 
tency tests of the current cosmological paradigm, which 
will continue with the E-ELT. A roadmap for these tests 
can be found in [49] . 

ACKNOWLEDGMENTS 

This work was done in the context of the project 
PTDC/FIS/111725/2009 from FCT (Portugal). C.J.M. 
is also supported by an FCT Research Professorship, con¬ 
tract reference IF/00064/2012, funded by FCT/MCTES 
(Portugal) and POPH/FSE (EC). 


[1] C. J. A. P. Martins, Phil.Trans.Roy.Soc.Lond. A360, 
2681 (2002), arXiv:astro-ph/0205504 [astro-ph] 

[2] J.-P. Uzan, Living Rev.Rel. 14, 2 (2011), 

arXiv:1009.5514 [astro-ph.CO] 

[3] T. Rosenband, D. B. Hume, P. O. Schmidt, C. W. Chou, 
A. Brusch, L. Lorini, W. H. Oskay, R. E. Drullinger, 
T. M. Fortier, J. E. Stalnaker, S. A. Diddams, W. C. 
Swann, N. R. Newbury, W. M. Itano, D. J. Wineland, 
and J. C. Bergquist, Science 319, 1808 (2008) 

[4] L. Amendola, A. C. O. Leite, C. J. A. P. Mar¬ 
tins, N. J. Nunes, P. O. J. Pedrosa, and 

A. Seganti, Phys.Rev. D86, 063515 (2012) 

arXiv:1109.6793 [astro-ph.CO] 

[5] J. K. Webb, J. A. King, M. T. Murphy, V. V. Flambaum, 
R. F. Carswell, et al., Phys.Rev.Lett. 107, 191101 (2011), 
arXiv: 1008.3907 [astro-ph.CO] 

[6] P. Molaro, M. Centurion, J. Whitmore, T. Evans, 
M. Murphy, et al, A. & A. 555, A68 (2013) 
arXiv:1305.1884 [astro-ph.CO] 

[7] H. Rahmani, M. Wendt, R. Srianand, P. Noterdaeme, 
P. Petitjean, P. Molaro, J. B. Whitmore, M. T. Mur¬ 
phy, M. Centurion, H. Fathivavsari, S. D’Odorico, T. M. 
Evans, S. A. Levshakov, S. Lopez, C. J. A. P. Martins, 
D. Reimers, and G. Vladilo, MNRAS 435, 861 (2013) 
arXiv:1307.5864 [astro-ph.CO] 

[8] T. M. Evans, M. T. Murphy, J. B. Whitmore, T. Misawa, 
M. Centurion, S. D’Odorico, S. Lopez, C. J. A. P. Mar¬ 
tins, P. Molaro, P. Petitjean, H. Rahmani, R. Srianand, 
and M. Wendt, M.N.R.A.S. 445, 128 (2014) 

[9] M. C. Ferreira, O. Frigola, C. J. A. P. 

Martins, A. M. R. V. L. Monteiro, and 

J. Sola, Phys.Rev. D89, 083011 (2014) 

arXiv:1405.0299 [astro-ph.CO] 

[10] A. Coc, N. J. Nunes, K. A. Olive, J.-P. Uzan, 

and E. Vangioni, Phys.Rev. D76, 023511 (2007) 
arXiv:astro-ph/0610733 [astro-ph] 

[11] A. Songaila and L. Cowie, Astrophys.J. 793, 103 (2014) 
arXiv:1406.3628 [astro-ph.CO] 

[12] D. Albornoz Vasquez, H. Rahmani, P. Noter¬ 

daeme, P. Petitjean, R. Srianand, and C. Ledoux, 
A. & A. 562, A88 (2014) 

[13] J. Bagdonaite, W. Ubachs, M. Murphy, and 

J. Whitmore, Phys.Rev.Lett. 114, 071301 (2015) 

arXiv: 1501.05533 [astro-ph.CO] 

[14] N. Kanekar, W. Ubachs, K. Menten, 


J. Bagdonaite, A. Brunthaler, et al, 
Mon.Not.Roy.Astron.Soc. 448, L104 (2015) 
arXiv: 1412.7757 [astro-ph.CO] 

[15] B. A. Campbell and K. A. 

Olive, Phys.Lett. B345, 429 (1995) 

arXiv:hep-ph/9411272 [hep-pli] 

[16] F. Luo, K. A. Olive, and J.- 

P. Uzan, Phys.Rev. D84, 096004 (2011) 

arXiv:1107.4154 [hep-ph] 

[17] M. C. Ferreira, M. D. Juliao, C. J. A. P. Martins, and 
A. M. R. V. L. Monteiro, Phys.Rev. D86, 125025 (2012) 
arXiv:1212.4164 [hep-ph] 

[18] A. Avgoustidis, G. Luzzi, C. J. A. P. Martins, and 


[19] 


A. M. R. V. L. Monteiro, 

JCAP 1202, 013 (2012) 

arXiv:1112.1862 [astro-ph.CO] 



V. L. Monteiro, P. E. Vielzeuf, and G. Luzzi, 


JCAP 1406, 062 (2014), arXiv: 1305.7031 [astro-ph.CO] 

[20] N. Kanekar, J. N. Chengalur, and 

T. Ghosh, Ap.J.Lett. 716, L23 (2010) 

arXiv: 1004.5383 [astro-ph.CO] 

[21] J. Darling, Astrophys.J. 612, 58 (2004) 

arXiv:astro-ph/0405240 [astro-ph] 

[22] M. Murphy, J. Webb, V. Flambaum, 
M. Drinkwater, F. Combes, et al ., 
Mon.Not.Roy.Astron.Soc. 327, 1244 (2001) 

arXiv:astro-ph/0101519 [astro-ph] 

[23] N. Kanekar, G. I. Langston, J. T. Stocke, C. L. Car- 
illi, and K. M. Menten, Ap.J.Lett. 746, L16 (2012) 
arXiv:1201.3372 [astro-ph.CO] 

[24] H. Rahmani, R. Srianand, N. Gupta, P. Petitjean, P. No¬ 
terdaeme, and D. A. Vasquez, MNRAS 425, 556 (2012) 
arXiv: 1206.2653 [astro-ph.CO] 

[25] S. Curran, A. Tanna, F. Koch, J. Berengut, 

J. Webb, et al, Astron.Astrophys. 533, A55 (2011) 
arXiv:1108.0976 [astro-ph.CO] 

[26] A. Weiss, F. Walter, D. Downes, C. Carrili, 
C. Henkel, et al, Astrophys.J. 753, 102 (2012) 
arXiv: 1204.5614 [astro-ph.CO] 

[27] R. Srianand, N. Gupta, P. Petitjean, P. Noter¬ 

daeme, and C. Ledoux, MNRAS 405, 1888 (2010) 
arXiv: 1002.4620 [astro-ph.CO] 

[28] L. Lentati, C. Carilli, P. Alexander, R. Maiolino, 

R. Wang, P. Cox, D. Downes, R. McMahon, 

K. M. Menten, R. Neri, D. Riechers, J. Wagg, 




F. Walter, and A. Wolfe, MNRAS , 721 (2013) 

arXiv:1211.3316 [astro-ph.CO] [39] 

[29] S. A. Levshakov, F. Combes, F. Boone, I. I. Agafonova, 

arXiv:1104.3361 [astro-ph.CO] 

J. Bagdonaite, M. Dapra, P. Jansen, H. L. Bethlem, 
W. Ubachs, et al, Phys.Rev.Lett. Ill, 231101 (2013) 

D. Reimers, and M. G. Kozlov, A. & A. 540, L9 (2012) 

arXiv:1311.3438 [astro-ph.CO] 

arXiv: 1203.3649 [astro-ph.CO] [40] 

[30] M. C. Ferreira, M. D. Juliao, C. J. A. P. Martins, and 

F. van Weerdenburg, M. Murphy, A. Malec, L. Kaper, 
and W. Ubachs, Phys.Rev.Lett. 106, 180802 (2011) 

A. M. R. V. L. Monteiro, Phys.Lett. B724, 1 (2013) 

arXiv:1104.2969 [astro-ph.CO] 

arXiv:1305.7515 [astro-ph.CO] [41] 

[31] P. Molaro, D. Reimers, I. I. Aga¬ 
fonova, and S. A. Levshakov, 

A. L. Malec, R. Bulling, M. T. Murphy, 
N. Milutinovic, S. L. Ellison, J. X. Prochaska, 
L. Kaper, J. Tumlinson, R. F. Carswell, 

European Physical Journal Special Topics 163, 173 (2008), 

and W. Ubachs, MNRAS 403, 1541 (2010) 

arXiv:0712.4380 

[32] H. Chand, R. Srianand, P. Petitjean, B. Aracil, [42] 

arXiv:1001.4078 [astro-ph.CO] 

J. Bagdonaite, M. T. Murphy, L. Kaper, 

R. Quast, and D. Reimers, A. & A. 451, 45 (2006) 

and W. Ubachs, MNRAS 421, 419 (2012) 

arXiv:astro-ph/0601194 

[33] I. I. Agafonova, P. Molaro, S. A. Levshakov, [43] 

arXiv:1112.0428 [astro-ph.CO] 

J. A. King, J. K. Webb, M. T. Murphy, and 

and J. L. Hou, A. & A. 529, A28 (2011) 

R. F. Carswell, Phys. Rev. Lett. 101, 251304 (2008) 

arXiv:1102.2967 [astro-ph.CO] 

[34] M. T. Murphy, V. V. Flambaum, S. Muller, [44] 

arXiv:0807.4366 [astro-ph] 

J. A. King, M. T. Murphy, W. Ubachs, 

and C. Henkel, Science 320, 1611 (2008) 

and J. K. Webb, MNRAS 417, 3010 (2011) 

arXiv:0806.3081 [astro-ph] 

[35] N. Kanekar, Astrophys.J. 728, L12 (2011), [45] 

arXiv:1106.5786 [astro-ph.CO] 

M. Wendt and D. Reimers, 

arXiv:1101.4029 [astro-ph.CO] 

[36] C. Henkel, K. M. Menten, M. T. Murphy, N. Jethava, 

V. V. Flambaum, J. A. Braatz, S. Muller, [46] 

Eur.Phys.J.ST 163, 197 (2008), 
arXiv:0802.1160 [astro-ph] 

M. Kozlov and S. Levshakov, 

J. Ott, and R. Q. Mao, A. & A. 500, 725 (2009) 

Annalen Phys. 525, 452 (2013), 

arXiv:0904.3081 [astro-ph.CO] 

[37] V. V. Ilyushin, P. Jansen, M. G. Kozlov, [47] 

S. A. Levshakov, I. Kleiner, W. Ubachs, and 

arXiv: 1304.4510 [physics.atom-ph] 

J. Bagdonaite, W. Ubachs, M. T. Murphy, and 
J. B. Whitmore, Astrophys.J. 782, 10 (2014) 

H. L. Bethlem, Phys.Rev. A85, 032505 (2012) 

arXiv:1308.1330 [astro-ph.CO] 

arXiv: 1201.2090 [chem-ph] [48] 

J. B. Whitmore and M. T. Murphy, 

[38] S. Muller, A. Beelen, M. Guelin, S. Aalto, J. H. 

Black, F. Combes, S. J. Curran, P. Theule, 

Mon.Not.Roy.Astron.Soc. 447, 446 (2015) 
arXiv: 1409.4467 [astro-ph.IM] 

and S. N. Longmore, A. & A. 535, A103 (2011), [49] 

C. J. A. P. Martins, Gen.Rel.Grav. 47, 1843 (2015) 


arXiv: 1412.0108 [astro-ph.CO] 




